function y = computeKinematicHardeningModulus(element,alpha)
H = element.KinematicHardeningModulus;
y = H*alpha;
end
